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, We present a new calculation of the tt-J/i/j dissociation cross sections within the Constituent 

^ ' Quark-Meson Model recently introduced. To discuss the absorption of J/il) in heavy-ion collisions, 

D ' we assume the J/ip to be produced inside a thermalized pion gas, as discussed by Bjorken, and 

introduce the corrections due to absorption by nuclear matter as well. We fit the absorption length 
' of the J/tj) to the data obtained at the CERN SPS by the NA50 Collaboration for Pb-Pb collisions. 

' Collisions of lower centrality allow us to determine the temperature and the energy density of the pion 

gas. For both these quantities we find values close to those indicated by lattice gauge calculations 
■ for the transition to a quark-gluon plasma. A simple extrapolation to more central collisions, which 

^ ' takes into account the increase of the energy deposited due to the increased nucleon flux, fails to 

, reproduce the break in J/t/j absorption indicated by NA50, thus lending support to the idea that 

. an unconfined quark-gluon phase may have been produced. This conclusion could be sharpened 

fS) ' by analysing in a similar way, as a function of centrality, other observables such as strange particle 

fSJ ' production. 

o : 

■<sj- . PACS numbers: 25.75.-q, 12.39.-x 

O ; 

p I. INTRODUCTION 

(D : 

\ In a recent paper Q the couplings of J /ip to D^*^ mesons (where by D^*^ we mean D or D*) and to tt's have been 

J> ■ computed in the framework of a Constituent Quark-Meson Model introduced in . The interest of this study stems 

' from the possibility that J /ip absorption processes of the type: 

^ ^ + J/V' -> i^W + ^ (1) 

play an important role in the relativistic heavy- ion scattering. Since a decrease in the J/ip production in these 
processes might signal the formation of Quark-Gluon Plasma (QGP), it is useful to have reliable estimates of the 
cross sections for the processes which provide an alternative way to reduce the J/ip production rate. Previous 
studies of these effects can be found in ■ In this paper we will address the problem of understanding how far we 
can go in explaining J /ip suppression data via an hadronic mechanism such as . Does the hot gas of pions formed 
after the heavy-ion collision provide a source of attenuation of J /ip antagonist to the standard QGP suppression? 

The relevant couplings needed to compute the cross section of process ^ are shown in Fig. ^ To compute the 
amplitudes besides the DD*tt couphng, Fig^, whose strength go'D-n bas been theoretically estimated and 
experimentally investigated 3, one would need also the JD^*'^D^*\ Fig^, and the JD^*^D'^*^7r couphngs, FigOJ;. 
In an effective Lagrangian approach the latter couplings provide direct four-body amplitudes, while the former enter 
the amplitude via tree diagrams with an exchange of a charmed particle D^*^ in the t-channel. 
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FIG. 1: Couplings involved in the tree— level calculations of processes ITl . 



These couplings have also been estimated by different methods, all presenting, in our opinion, different shortcomings. 
For example the use of the SU (4) symmetry puts on the same footing the heavy quark c and the light quarks, which 
is at odds with the results obtained within the Heavy Quark Effective Theory, where the opposite approximation 
uic S> ^QCD is used. Similarly, the rather common approach based on Vector Meson Dominance should be considered 
critically, given the large extrapolation from = to = 'm^j/^ that is involved. A different evaluation, based on 
QCD Sum Rules can be found in and presents the typical theoretical uncertainties of this method. It is useful 
to have diverse, although inevitably model-dependent, calculations to assess a reliable range of theoretical results. A 
comparison of our results with the cross sections obtained by other methods is given below, so to provide a further 
indication of the theoretical uncertainties (the theoretical uncertainties intrinsic to the model used in this paper have 
been discussed in U). 

To apply our results to J /ip production in heavy-ion collisions, we follow the picture introduced by Bjorken some 
time ago |7|. Seen in the centre-of-mass frame, the colliding nuclei appear as Lorentz-contracted pancakes which 
traverse one another in a short time t (say, 2 fm/c). After crossing, the nuclei leave behind them a fireball with 
essentially vanishing baryon number, where the produced J /ip is sitting (particles in the fireball are often referred to 
as comovers). 

For low-centrality collisions, the fireball should be, to a good approximation, a pion gas that quickly goes to thermal 
equilibrium with a temperature corresponding to the energy deposited by the colliding nucleons in the central region. 
The energy density estimated by Bjorken is in the range of 0.3-3 GeV/fm'^ and is proportional to the number of 
nucleons per unit area that participate in the collision. 

Our cross section calculation allows us to compute the absorption length of the J/ip in the pion gas, which we 
tentatively approximate with a perfect gas with vanishing chemical potential (at high temperature, say T > 150 MeV, 
there are enough inelastic collisions in the pion gas to provide chemical equilibrium as well). 

To obtain a more accurate value of T, we introduce also the corrections due to J/ip absorption by nuclear matter, 
using the absorption cross sections obtained recently from p — A collisions mL 

We have fitted the absorption length to the data obtained by the NA50 p\ collaboration at the SPS, which show 
a clear exponential behaviour as a function of the linear dimension of the collision region, Z, up to / ~ 5 fm. We find 
not unreasonable values for the temperature (T = 225 McV) and for the energy density (e ~ 0.32 GeV/fm"^). The 
absorption length decreases very steeply with increasing temperature. 

The energy density is perhaps on the low side with respect to a priori expectations based on the Bjorken formula. 
On the other hand, the temperature is somewhat higher than the critical temperature computed in lattice calculations 
(see for a recent review). The difference, if any, could be due to either an underestimate of the cross section, or 
to the perfect gas approximation being too crude, or to both. Large departures from perfect gas behaviour are 
found in numerical simulations below the critical temperature P). Further work to improve on the perfect Bose 
gas approximation and obtain a better calibration of the crucial parameters of the fireball, temperature and energy 
density, is needed. 

To extrapolate to higher centrality, we keep into account the increase of the average nucleon number per unit area 
for decreasing impact parameter. This, in turn, makes the temperature and the deposited energy increase while the 
absorption length decreases, leading to a downward bending of the J/ip production as a function of /. The behaviour 
we find, however, is too smooth to reproduce the rather sharp break shown by the NA50 data for / larger than 5 fm. 
This suggests that the simple pion gas description ceases to be valid at these values of I, thus lending support to 
the formation of a new phase. Our analysis implies that it would be useful to analyse in the same way, i.e. as a 
function of I, the data obtained for other observables such as the production of strange particles. Given the rather 
low temperatures that we find in the low-centrality region, we would expect very few strange particles to be produced 
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in the central region as well. This would not apply to the unconfined phase, where the Boltzmann suppression of 
strange quarks with — 150 MeV would not be operative for temperatures around 225 MeV. 



II. CROSS SECTIONS 

The calculation of the total cross sections a{7rJ/^ £)(*)£)(*)) proceeds through the evaluation of the amplitudes 
listed in the Appendix and the two-body phase-space integration: 



where s is the Mandelstam variable s = (p^ R is the phase-space interval R = (p^j'^'P™'*'^)- overall 

factor of 2 is included to count all isospin quantum numbers in the final state. The laboratory frame is considered 
with high- momentum massless pions colliding on J/ip at rest. En energy threshold of about E^^ — 800 MeV is required 
to open the reaction channel. 

For example, in the case of the DD final state, using the notation introduced in 0, we have A = A^a + ^2a + ^3a, 
where: 

Ala = I JZTJ^ e('7,P4,P3 -P2,a)n [ps ~ P2)[P2ju 

Alb = « 7j2 e(^,P3,Pi -P3,a)n [p3-Pl)[P2)y 

Ale = i—r-r — e('7,P2,P4,P3), 

where rj is the polarization vector of J ftp, n^''(<z) is the sum-over-polarization-tensor 

^^''(9) = ^-5^^ (2) 

M being the mass of the vector particle. The notation e(/i, z^, p, g) = tf^upaP'^'f is used. The coupling constants 
g-nDD' , gj /■4)DD* T gjDDi^ have been discussed in Q (the notation gjoDrr = go is used), and are: 
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g-rrDD' — 



-1\ 



gjH'DD' = 4.05 (GeV 
gjDD^ = , , q^-p,,^ (GeV ). 
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where Ai = 1.1 GeV and A2 = 0.8 GeV. Attempts to compute quadrilinear couplings using SU{4) symmetry can be 
found in The dependence of the cross sections on -/s for the DD and DD* final states are shown in Fig.|21 

The amplitudes for DD*,D*D* final states are given in the Appendix. The results we obtain can be compared 
with those shown in Fig. |21 extracted by the authors of Note the sharp rise of DD* and D*D*, which is due to 
5-wave production. 



III. THERMAL AVERAGE 



In a first instance we neglect J/tp absorption by the nuclear matter, and focus on the effect of the fireball left over 
after the nuclei have passed each other (i.e. the comoving particles in the frame of the target nucleus). The quantity 
we are interested in is the J /ip suppression as a function of the linear size I of the fireball. A J /ij) produced in any 
internal point of a fireball of diameter I has to travel on average a distance ~ 6/lOZ to reach the boundaries of the 
ball. On the other hand the mean free path of the J/i/' is A ~ l/ptr, so the attenuation can be defined as: 



A{x) = exp [-x{p ■ cr^j/^^DMDM )t\ , 



(3) 
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FIG. 2: Cross sections c^j/^ — » {DD, DD* , D* D") as functions of We show the error band for the DD final state resulting from 
the theoretical uncertainty in the determination of the couplings yjj. The error bands on DD* and D* D* are similar. 
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FIG. 3: Cross sections computed with various approaches: QCD sum rules (band), short-distance QCD (dotted line), meson-exchange 
models (dot-dashed line), non-relativistic constituent quark model (dashed line) .10|. 



with X ~ 6/lOZ. 

As a simple approximation, the thermal average (...)t is taken in a perfect Bose gas of charged and neutral pions 
with vanishing chemical potential: 



3 f°" 

{P ■ (^■nJ/i,^D<.*)D<.*))T = TT^ f ^E^ 



For reference, we give also the number density of the relevant pions at temperature T 



(4) 
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El 



(5) 
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where i?* is the threshold energy required to open the reaction channel, and the total energy density: 
We report in Table I the values of p(T), e(T) and (p • (t)t for different values of the temperature. 



r (MeV) 


P (fm-^) 


e(T) (MeV/fm'') 


(P(j)t (fm) ^ 


205 


0.0467 


215 


0.037 


215 


0.0629 


261 


0.051 


225 


0.0828 


313 


0.070 


235 


0.11 


373 


0.093 


245 


0.13 


441 


0.12 


255 


0.17 


518 


0.16 


265 


0.20 


605 


0.20 


275 


0.25 


702 


0.25 



TABLE L Values of p{T), e{T) and {p ■ (j)t for different values of the temperature. 




FIG. 4: The overlap region of two colhding equal nuclei. Dimensions I and L determine the transverse size of the pion fireball and the 
longitudinal size of the nucleon column to be traversed by the J/^j), respectively. 



To connect to experimental data, we observe that, in the approximation of sharp-edged nuclei, the linear size of 
the fireball, I, is given by the transverse dimension of the region where nuclei overlap (see Fig. which in turn is 
related to the impact parameter b according to: 

l = 2R-b . (7) 
The NA50 Collaboration has studied the attenuation of J/ip by studying the ratio: 



^'^''P* " a(Drell - Yan) " ^ ' 

as a function of the impact parameter Q . Using Eq. {Tj), we can fit our Eq. © to the NA50 data, plotted as 
a function of I. The result is shown in Fig. |S1 for different values of the fireball temperature. A good fit can be 
achieved in the region of lower centrality, up to Z ~ 5 fm, thereby determining the fireball's temperature. Observe 
that peripheral interactions agree with extrapolation from p- or d- induced reactions (the points at Z ~ 0). From Fig.[Sl 
we find T ~ 255 MeV. 

To be more precise on the temperature, however, we must consider the effect of nuclear absorption, due to the 
nucleon column density, which the J/V' has to traverse during the initial phase of the collision. Neglecting the J/V' 
transverse momentum, the nuclear absorption factor is: 



-^nucl.abs. — CXp [ L/ • Pnucl ' f^nucl] 



(9) 
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L is approximately equal to the longitudinal size of the overlap region of the colliding nuclei (Fig. 0J. To avoid 
unwanted distortions, we have interpolated the values of L given by NA50 as functions of the impact parameter b to 
obtain L as a function of our length I via Eq. Q. A determination of the nuclear absorption cross section has been 
obtained by NA50 from the study of J /ip production in p-A collisions. They find : 

cr„„^i = 4.3 ± 0.6 mb (10) 

This cross section is smaller than what would explain the data in the low-centrality region with nuclear absorption 
only (this requires (Tnuci — 6 mb 0) and points to the importance of dissociation induced by the comoving particles, 
which in fact turns out to dominate. In conclusion, we compare the experimental quantity Eq. ^ to the expression: 



Ali(0 = • exp [-L(0 ■ Pnucl ■ O-nucl] • exp [-6/10 • • cr)T] , (11) 

where C is a normalization factor; L and I are obtained from the impact parameter 6 following the NA50 prescription 
and Eq. {Tj), respectively. We use pnuci — 0.17 fm"^, anuci — 4.3 mb. 

The results are shown in Fig. |H1 The fireball temperature that gives the best fit to the low-centrality data, up to 
I ~ 5 fm, is now T ~ 225 MeV, and the corresponding energy density is e ~ 0.32 GeV • fm^'^. We have studied how the 
theoretical errors in the coupling constants determined in u!| aff'ect the determination of T. We find T = 225±15 MeV. 

Going towards higher centrality, we expect an increase in the energy density of the fireball. One effect that is easy 
to take into account is the increase in the surface density of the nucleons participating in the collision. The energy 
density of the fireball is, in fact, proportional to the factor 0: 

^^|...(V.) (12) 



3 arccos(6/2i?)-(V2i?)v/l-^V4i?' 

with R = A^/'^ro, the nuclcon radius ro = 1.1 fm [Tll |. 

Starting from — 5 fm and Tq = 225 MeV, the geometrical factor g{2 ~ l/R) is used to extrapolate the variation 
of temperature with I according to: 



rm-T f 9i'2-i/R) \ 

^^'^"^'[9i2-lo/R)) 



1/4 

(14) 



Once we know T{1) we compute the corresponding attenuation function defined as: 

J^(/) = C • exp [-i(/)/9nuci • CTnuci] exp [-6/10 • r (p • cr)T(o] ■ (15) 

Fig El shows the final best prediction of the conventional absorption effects due to the pion gas fireball and to the 
nuclear matter (solid curve). 

Fitting the data of low-centrality, we are able to determine with some accuracy the temperature of the fireball, 
assumed to be an ideal pion gas. The temperature is remarkably close to, perhaps somewhat higher than, the critical 
temperature computed in lattice calculations 0, '^-^ 180 MeV. We recall that in SU(3) gauge theory, Tc = 260 MeV 
while lattice results with Uf — 2 suggest Tc = 170 MeV. If the same high temperature limit is assumed for > 
and SU(3), then for n/ = 2 or 3, the ip' breaks at 190-200 MeV and the J/i/' itself would not break up until T > 2Tc 
[l_^. Departures of the estimated temperature from the real one could be due either to an underestimation of the 
cross section, or to the inadequacy of the perfect gas approximation, or to both. Large deviations from the perfect 
gas are found in numerical simulations below the critical temperature Q. Further work to improve on the perfect 
Bose gas approximation and to obtain a better calibration of the crucial parameters of the fireball, temperature and 
energy density, is needed. 

The energy density can be compared with the formula given in ref . : 

A{b) (dE\l 

A{b) and S{b) are the number of nucleons that participate in the collision and the overlap area, respectively, as 
functions of the impact parameter b; dE/dy is the energy deposited per unit rapidity and ct, the longitudinal length 
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of the fireball, is related to the time t it takes for the nuclei to separate completely. The energy deposited per unit 
rapidity is estimated as follows : 



dy dy \ N^n 
For central Pb-Pb collisions, one has: 



dE dN,h ^ ]{E)^?,y. 1.5 x 400 MeV = 1.8 GeV 



(17) 



For Pb and L = 5 fm we find g — 0.7. In conclusion, for Pb-Pb collisions we find: 



(18) 



e = 0.95 GeV/fm^ ["^ 



e 1.35 GcV/fm''' 



(L = 5 fm) 



2 fm 

ct 



{L = 2R= 13 fm) 



(19) 
(20) 
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FIG. 5: The attenuation function A given in Eq. @ as a function the linear size of the fireball. Lines correspond to a 
fixed temperature of the tt gas expressed in MeV (isotherms). Experimental data: diamonds^ p-{p,d) (NA51), boxes=Pb- 
Pb (NA50). The value of I for each experimental point is obtained from the measured value of the impact parameter b via 
Eq. ||7|l. Lower-centrality interactions agree with extrapolation from p- induced reactions. 



IV. CONCLUSIONS AND OUTLOOK 



Our final results are displayed in Figs. El and [3 In the lower-centrality region, where we expect to be still in the 
deconfined phase, we find T = 225 ± 15 MeV corresponding to the energy density e ~ 0.32 GeV/fm^. The value of 
the energy density is not inconsistent with the Bjorken formula. The temperature, however, is somewhat higher than 
the critical temperature computed in lattice calculations. 

The extrapolation to higher centrality keeps into account only the increase of the average nuclcon number per unit 
area (an admittedly crude approximation). This leads to a downward bending of the J/^ production as a function of 
I. However, the behaviour is too smooth to reproduce the rather sharp break shown by the NA50 data. 

Our analysis suggests that the simple pion gas description ceases to be valid at higher values of I, thus lending 
some support to the formation of a new phase. To confirm or disprove this, it would be crucial to analyse in the same 
way, i.e. as a function of I, the data obtained at the SPS and RHIC for other observables, such as strange particles 
production. In a pion gas, given the relatively low temperatures that we find in the low-centrality region, we would 
expect very few strange particles to be produced in the central region. This would not apply to the deconfined phase, 
where the Boltzmann suppression of strange quarks with nig = 150 MeV would not be operative for temperatures 
around 225 MeV. 
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FIG. 6: The attenuation function A given in Eq. 1111 as a function of the linear size of the fireball, computed with Pnucl = 0.17 fm~^ 
and (T,yyr.i = 4.3 mb. Lines correspond to a fixed temperatur e of the tt gas expressed in MeV. Experimental data; Pb-Pb collisions from 
NA50 13. Diamonds represent p — (p,d) (NA51) data points Il3l . 
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FIG. 7: Complete calculation of the attenuation function. T = 225 MeV, Pnucl = 0.17 fm~^ and (Tnud = 4.3 mb. The solid line includes 
the effect of the increase of nucleon density per unit surface with decreasing impact parameter. The dashed line, given for reference, is the 
same as in Fig. 6. Data from NA50 Pb-Pb collisions j^. 



Appendix 

We report the amplitudes corresponding to the tree-level diagrams for the DD* and D*D* final states. The 
diagrams and the corresponding amplitudes are as follows (dashed lines for spinless particles, continuous lines for 
spin-1 particles); we use the same notation as in 0. 




J/iP D J/i, D J/tP D 

(la) (lb) (Ic) 
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In the case in which the final state is DD* : 



= t-Ml, (P2 • eb3))(p4 • e(pi)) 

e(e(P3),P2,P3, 1^)11'''' (p3 -P2) 

((Ps - Pi) ■ e(Pi))e,i(P3) + • e(p3))p3^)x 

n^^bs - pi, Ml,.) + ((P3 - Pi) • e(p3))e.(pi)x 

n-''(p3~Pi,M|,.)) 

-42d = -P2e^(Pl) (.93M^^(P3i'P4^(P4 • e(P3)) -P3mP4i.(P4 • e(P3)))- 

92{piv^ix{p3) + Pzv^^x{pz)) + 9i{9iiv{pi ■ e(P3))+ 
2p3,*e!.(P3))) 



The amplitudes for the D*D* final state are: 



^3a = j ^^^^^^^^^^ e(e(ffi), e(p4),P3 - P2,Pi){P2 ■ e{pz)) 

^35 = ^ ^'"^ ^ •'^,'^^2^^* g(g(Pl)^P3-P2,e(p3),P3)(p2-e(p4)) 
M — M£) 

. 9-nD'D'9j/rpD'D* , , , 

^ t-M^ {e{P3 -P2,r,p3,e(P3)) X 

n^"(p3 -P2,MD.)(P4A(ebi) • e(p4)) - 
eA(pi)((p3 - P2) • e(p4)) - eA(p4)(p4 • e{pi)) + 
eA(p4)((p3 -P2) • e(pi))) 

Asd = 772 e(P4,e(P4),P3 -Pi,P) X 

n^^(P3 - Pi,MD' )i-e{p3) ■ e(Pi)P3A + (P3 • e(Pi))eA(P3) + 

{(P3 -Pi) • e(Pi))eA(P3) - 
((P3 -Pi) • e(P3))eA(pi)) 

Ase = iMDP2u^,j.(pi) ^-^e(/x,e(p3),e(p4),P3)P3 - ^e(/x,e(p3),e(p4),P3)P4 + 
^i-^^e(/i,p4,e(P4),P3)e''(P3) + -^e(/i, e(p3),P4, e(P4))P3 + 
-^e(j/,p4,e(p3),e(P4))P4 + ^e(M,P4, e(P3), e(P4))P4 + 
^^^2^ e{v,pi, e(p3),p3)e^(p4) - 9Ae{lJ,,v, e{p3),e{pi))- 

^^^^^r^e(p4,e(P3),e(P4),P3)P3P3 + ;^e(P4, e(P3), e(p4),P3)P4P4- 

-^(2e(z^,e(p3),e(p4),P3)P3 - e(P4,z^, e(P4),P3)e''(P3)+ 
2e(/i,p4,e(p3),P3)e''(P4) - 5'"'e(P4,e(P3),e(P4),P3))) 
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